Ground-state factorization and quantum phase transition in dimerized spin chains 
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We study the occurrence of ground-state factorization in dimerized XY spin chains in a transverse 
field. Together with the usual ferromagnetic and antiferromagnetic regimes, a third case emerges, 
with no analogous in translationally-invariant systems, consisting of an antiferromagnetic Neel-type 
ground state where pairs of spins represent the unitary cell. Then, we calculate the exact solution 
of the model and show that the factorizing field represent an accidental degeneracy point of the 
Hamiltonian. Finally, we extend the study of the existence of ground-state factorization to a more 
general class of models. 



PACS numbers: 03.67.Mn, 75.10.Jm, 64.70.Tg 



The study of zero-temperature critical phenomena in 
quantum magnetic systems represents since long time 
a major research subject lii^ii^ii^ In particular, the XY 
model is a very rich source of information about the 
quantum behavior of spin chains because of the avail- 
ability of an exact analytical solution. '^■^ During last 
years, the main interest about spin chains concerned 
the relationship between quantum phase transitions and 
entanglemen1j^iii2ii^. Among a number of interesting 
properties of such systems, it is worth citing the ex- 
istence of special values of the external magnetic field, 
the parameter which drives the phase transition, which 
give rise to ground state factorization, discovered by Kur- 
mann et al. M- This phenomenon has been observed in 
two-dimensional lattices through quantum Monte Carlo 
methods''''' and fully analyzed by Giampaolo et al}^ in a 
recent publication , where the factorizing field has been 
determined for a quite general class of models. They de- 
veloped an appropriate measure of entanglement which 
vanishes at the factorizing point. So far, the existence 
of a factorized ground state has been predicted only in 
translationally invariant Hamiltonian models. 

Moreover, critical properties of physical systems are 
discussed by taking the thermodynamic limit from the 
beginning. On the other hand, the knowledge of a finite- 
size solution clarifies important aspects of this limit. For 
example, it is known that the quantum phase diagram of 
the XY chain in a transverse field exhibits two different 
symmetry-broken regions characterized by different be- 
haviors of two-body correlation functions. This dissim- 
ilarity has a microscopic origin easily understood in the 
finite-size case. Besides these considerations, the study 
of finite systems is relevant by itself for the realization of 
mesoscopic qubits of contemporary interest. 

In this Rapid Communication we discuss a finite- 
size dimerized XY spin chain in a transverse field, and 
analyze ground-state properties. The interest about 
such system, belonging to a more general class of 
modelsf^ii^ii^ is motivated by experimental work on qua- 
sicrystals and quasiperiodic superlattices--iS*i^ First of all, 
we show that the model admits the existence of a factor- 
ized ground state, and then discuss the exact solution. 
The factorizing field turns out to be an accidental de- 



generacy point of the Hamiltonian and falls on a border 
surface between two regions that, in the thermodynamic 
limit, are characterized by different symmetry-breaking 
mechanisms. Furthermore, we will be able to detect the 
conditions for the existence of ground-state factorization 
also in a more general class of dimerized chains which are 
the generalization of the model discussed in Ref J^. 

We start our discussion by considering a nearest- 
neighbor dimerized chain of an even number N of spin 
1/2: 
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with cr^ 



= af, being af the ath Pauli matrix (a — 
X, y, z). Without loss of generality, we will limit the anal- 
ysis to positive fields. The Hamiltonian of Eq. ((T|) can be 
recast in the form of a sum over two-body Hamiltonians 



W = Erii'(<+WrO, where 
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with hi and such that h = hi + h2. 

The central feature of Tij is the invariance under ro- 
tations of TT around the z axis. This is formalized by 
the vanishing of the commutator [T-l[^\Vi ^^] — 0, where 

T'l^'' = cr2i_2j^i<^2i-i+i the parity operator, since its 
eigenvalues are -1-1 or —1, according to the number of 
down spins in the z direction being even or odd. The 
above commutation relation then requires also the eigen- 
states of Tij*-* to have definite parity. The problem is 
to establish whether there exists a set of the Hamilto- 
nian parameters such that the ground state is of the 
form l^*) = Notice that, if this is the case, as re- 

marked in Ref.""^", \'4'2i-2+i)\i^2i-i+i) must be the ground 



state of The problem is then reduced to find the 

conditions under which Tij'"* admits a factorized ground 
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state. Now, since 'H.f' is not diagonal in the tr^ basis, if 
|V'2!-2+i)|'02/-i+i) has to be the ground state, each fac- 
tor must be of the form |?/;;) = cosi/'; IT) +sini/); ||) , with 
ipi ^ 0,7r/2. This, in turn, imphes that in the factor- 
ized ground state the parity symmetry is broken and, 
therefore, looking for the condition on the parameters of 
the Hamiltonian leading to factorization of the ground 
state amounts to looking for the condition leading to 
the degeneracy of the even and odd lowest lying eigen- 
states, without invoking the vanishing of entanglement 
indicatorsJ^ 



There are three different physical scenarios to be con- 
sidered: (i) both Ji and J2 are negative (ferromagnetic 
case); (ii) both Ji and J2 are positive (antiferromagnetic 
case) ; (iii) Ji > and J2 < or vice versa (hybrid case) . 
In any of these cases, factorization appears if and only 
if 71/Ji = 72 / J2 — i.e., only in the presence of per- 
fect dimerization of the longitudinal part of 7i, the fac- 
torized point falls in hp — [( Ji -I- J2) /2] ^1 — 7^, and 



tan?/;2/-i+i = ± [(l - VT 



1/2 



(i) In the ferromagnetic case, we find ilj2i-i+i — 
ip2i-2+i- As expected, the factorized state is fully aligned 
along two possible directions: \"^) = 

(ii) In the second (antiferromagnetic) case, we find 
'02i-i-(-i = —'4'21-2+i- Then, alternate directions for the 
spin determining a Neel-type ground state are observed: 



I*) = (S> 



21+ 



1) |V'^+2>- 



(iii) If both ferromagnetic and antiferromagnetic fac- 
torized ground states are of the same kind of those ob- 
tained in the homogeneous XY chaiufiS the third (hy- 
brid) case shows up an original character. Indeed, by 
assuming, for example, Ji < and J2 > we find the 
constrains ip2i-i = "^21 and ^l)2i = —'^21+1- As a conse- 
quence, the factorized ground state assumes the structure 

I*) = ®tt^~'\^ti+,)K+2)Wi+,)\i^l+,)- Thus, 
we obtain an antiferromagnetic Neel-type ground state, 
whose unitary cell is represented by a pair of spins. An 
additional requirement for the existence of the FP in this 
case is that iV/4 must be an integer number in order to 
avoid frustration effects. 



Now, we discuss the general exact solution of Hamil- 
tonian ^ and enlighten the role played by the factoriz- 
ing field. The diagonalization method is given in RefJ^. 
We discuss explicitly the finite-size limit J^ii^ The first 
step is the introduction of the Jordan- Wigner transfor- 
mation, mapping spins into spinless fermions,^ defined 

through CTf = 1 - 2c\cu = X\j<i (l - ^A'^i) ^i, and 
CT;^ = Ylj^i ( 1 — 2cjci ] cj, which leads toH — T-Cq—VHi, 
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where the parity operator isV — 11/1=1 ~ ^"^z ■ Since 
\H,V] = 0, all eigenstates of H. have definite parity, and 
we can proceed to a separate diagonalization of Ti in 
the two subspaces corresponding to V = ±1. Then, 
the complete set of eigenvectors of H will be given by 
the odd eigenstates of Ti.^ = Ho + Hi and the even 
eigenstates of H'^ — Ho ~ Hi. Both for and H~ 
the diagonalization can be performed through the divi- 
sion of the lattice in two sublattices: C2i-i = a;, and 
C21 — bi, and with the help of two separate Fourier 
transforms a; 



{N/2)-'/^j:^akexp[-t^] and k 

3xp[-i^], where fc = 
1 in W-, and k = 1/2, 3/2, .. . {N/2) -1/2 in H+ , getting 



(iV/2)-^/' Efc 6fcexp[-z^], where fc = 0, 1, . . . (iV/2) - 



k 

- /i^(2-24afc-26[.6fc) , (5) 



with Jfc = Ji + J2 exp[-z^], 7fe = 71 - 72 exp[-i^]. 
We remark that the difference between 7i+ and con- 
sists in the different set of allowed values of fc. Finally, we 
introduce a generalized Bogoliubov transformation con- 
necting Cfc, aLfc, bk, feLfc, and obtain two different kinds of 
quasiparticles, 

H^^Y. ^i"' {''^i^^ - 1) +E (2el.a - 1) , (6) 



where the eigenvalues, belonging to two branches sep- 
arated by an energy gap, are, for fc ^ 0, 7V/4, 



A' 



(±) 



\/rk ± V^, with Tfc = Ah'^ + \Jk\' 



and Sk = \Jk\'{l6h^ + 2\^k\') + Jhf + Ipt^ 
and An 



4/i2 + (71 - 72)' ± (Ji + J2), A** 



N/i 



AK^ + (71 +72) ± (Ji — J2). Let us assume, for ex- 
ample, (Ji + J2) > and Ji > J2. Then, each A^, is pos- 
itive, with the exception of Aq~', A^^^ (both of them are 

eigenvalues of H~), which can assume negative values, 

f 1) / 2 2 
respectively, for h < h}-,' — -W ( Ji — J2) — (71 +72) /2 

and for h < h^^'> = \j [Ji + J2f ~ {n - 12? /2. Let 
us assume also, without loss of generality, /i^^ < 



3 



Due to the foregoing considerations, the ground state of 
Ti"*" is its vacuum, and the corresponding eigenvalue is 

a1''\ As for H", the low- 



P+ _ _Y-(JV/2)-(l/2)^- 
^0 — Z^fc=l/2 Z^iy= 



est energy is £^q 



Ml 



while the 



ground-state structure depends on h. The ground state 
has two quasiparticles on modes fc = 0, N/4 for h < 

(2) 
C ' 



one quasiparticle on the mode fc = for h^^? < h < h 

,(2) 



and finally, it is the vacuum state for h > h)j . 

In order to identify the ground state of we must 
compare the lowest eigenvalue of belonging to an 
even eigenstate (iJ™*^") with the lowest eigenvalue of 
belonging to an odd eigenstate {Eq'^'^). Since the vac- 
uum state is even, we immediately state that lowest even 



eigenvalue of Ti is i?. 



even 




As far as the lowest 



odd eigenvalue of Ti is considered, only inside the re- 
gion where there is one quasiparticle (odd 



odd 



In fact, outside 



number of excitations), E^ 
this region, Eq belongs to even eigenstates (vacuum or 
two-quasiparticle state), and we must look to the first 



excited state of H. 
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and E^'^'^ 



Thus, T^°'^'^ 
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{h < h 



c) 



E,{h > 
E^[h<h^^'^)+K^^l. 



In the thermodynamic limit, the sum over k becomes 
an integral, and the vacuum energies of 7i+ and 7i~ are 



identical: Eq 



En 



Thus, for h^f^'^ < h< h')^^ 



rpodd 

^0 - 



and spontaneous symmetry breaking comes out 
Outside this range, the ground state has definite par 



TTieven 
^0 



ity, being the energy gap AE = 2Aq ^ tor h > h^' and 
AE — 2A^^I^ for h < /i^-*. Thus, although phase transi- 
tions take place only for macroscopic systems, the change 
in the energy sign of quasiparticles can be used as precur- 
sive property, allowing one to characterize the transition 
before performing the thermodynamic limit. 

As in the homogeneous chain, ^ in the presence of spon- 
taneous symmetry breaking, two-body correlation func- 
tions can decrease monotonically or oscillate as a func- 
tion of the spin distance depending on the Hamiltonian 
parameters. Recently, de Lima et ali^ found different 
regions in the phase diagram separated by hypersurfaces 
(collapsing into hues when 71 /Ji = 72 /J2)- We show 
that the first one of the separating lines corresponds to 
the factorizing field. Indeed, the different behavior of the 
correlation functions derives from different symmetry- 
breaking mechanisms. By analyzing lowest odd and even 
eigenvalues of TL in the symmetry broken region for any 
finite N (Fig. [T]), we observe a series of N /2 intersection 
point hi. The existence of such points has been discussed 
for the homogeneous chain in Refs<^ ,ancl2i and is re- 
sponsible for magnetization jumps of Refi^^. If 71/ Ji — 
12/ the first crossing point (for decreasing fields) is at 

h = hp for any N. In fact, Yl,k=o~^ i^p) = 

J:kilj2^'^:^=+^k^ (hp) = N{J, + J2). That is, h = 
/li?, is an accidental degeneracy point. If 71/Ji ^72 /J2, 
there are not fixed points. However, with the increasing 
of iV, all the crossing point are confined inside the re- 
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FIG. 1: The energy difference Eq — Eg between the two 
ground states of Ti.^ and Ti.~ . Inside the range {/i^-* , /i^' } 
these two states are the lowest eigenstates of Ti. The dimen- 
sionless Hamiltonian parameters are Ji = 1, J2 = 0.4, 71 = 
0.42,72 = 0.168. These choice is compatible with the exis- 
tence of the factorized ground state (hp = h* (2)), which is 
found to be around h — 0.63. The other relevant parameters 
are ~ 0.06, h* (1) ~ 0.27, and ~ 0.69. The light gray 
line refers to a chain of 20 spins, the gray line corresponds to 
24 spin, while for the black curve a chain of 28 spin has been 
used. In the inset we plot a detail of the bigger picture in the 
case = 24 (hght gray) and N — 28 (black) to enlighten the 
level crossing. As expected, we found N/2 intersection points 
for each curve. 



gion limited hy h* {!) = \J ( Ji — — (71 — 72)^/2 and 

h* (2) = ^ ( Ji -I- - (71 -I- 72)^2- These critical val- 
ues define the separating surfaces of Refiii. It is worth 
noting that when a factorized field exists, hp = h* (2). 



In the thermodynamic limit, this kind of structure im- 
plies two different symmetry breaking mechanisms. As 
one can deduce from the results of Fig. [U for h* {!) < 
h < h* (2), as ^ 00, the set {hi} of the degeneracy 
points becomes a denumerable infinity, i.e., an infinite 
number of crossing points appears and the two lower 

energies coincide, while for /i^^"* < h < h* (1) and for 

(2) 

h* (2) < h < h)j there is the usual symmetry breaking 
due to the vanishing of the gap, i.e. the ground state has 
a definite parity for any finite N, but this difference goes 
to zero with 1/A^. This is the microscopic mechanism 
responsible for dissimilar two-body correlation functions. 
Then, like in the homogeneous case^ as the factorizing 
field is reached, correlation functions change character. 



To conclude, we extend the study of ground-state fac- 
torization to a larger class of XYZ dimerized long-range 
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spin models: H = J^r'^r, with 



a—x,y 
N 



(Af/2)-l N/2 

Eja a a , \ ^ ja a a 

'JrS"2l-l"2l-l+r ' r ,2" 2l" 2l^r 

l = \ 1 = 1 



(7) 



i=\ 

where J"j are the dimerized couphng constants between 
spin pairs at odd distance r. The existence of alternate 
coupling on even distances cannot be univocally intro- 
duced. However, we could consider homogeneous cou- 
pling for such distances. Notice that, assuming J"]^ = 
J^2 foi' ^iiy recover the class of models considered 

by Giampaolo et al. in Ref. For the sake of clarity, 
we will consider explicitly the case of a one-dimensional 
lattice, but the generalization to higher dimensions is 
straightforward. 

Let us first note that, to circumvent frustration effects, 
(N/r) has to be an integer number for any value of r 
appearing in 7i. As already done for the short-range XY 
chain, we rewrite Tir as a sum of two-body Hamiltonians, 
with 



n 



l,r 



rA'^2l-2+i" 2l-2+i+r 



cr-y 



rA'^2l-2+i"2l-2+i 



cr-y 



TV 

rS" 2l-2+i" 2l-2+i+r 
^i,r (o'2/-2+i + ^2l-2+i+r) J 

(8) 



and hi^r, ^2,r such that h = (/ii,r + h2,r)- 

The simplest case to study is the full ferromag- 
netic picture, where all coupling constants are neg- 
ative. By following the same procedure introduced 
above, we proceed to calculate the ground-state ener- 
gies for the two parity subspaces in each of two kinds 
of dimers {'Hf]^ and and force symmetry break- 

ing. The factorized ground state of each of such dimers 

(cos4'^|r)±sinVi'^|i))^ ® 
, with 



amounts to be l^* 



Lr 



^cos#|T)±sin4*^|i) 



l+r 



2 - 



tan^ Vr 



jx _ TV 

r.t r.i 
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, The existence of a globally factorized ground state im- 

(i) 

plies that the angle ipf- must be the same for any of 



the dimers involving each site, i.e. it has to be in- 
dependent both on r and i. This result is achieved 
by the following conditions: = ^Jr,i^ ^^^^ J"j — 

7rJf j. The value of the factorizing field is then hp = 

{1 + k) sJ{Jl - J'f ) (J'f - J'f), where the J^" are the 
global interactions along different axes: J"" = ^ri- 

For the full antiferromagnetic case, let first assume 
the nearest-neighbor interaction parameters such that 
(jf j 4- Jfj) > 0. Then, we expect for the factorized 

ground state the structure l^*) = \'^2i+\l |V'2;+2)- 

The range-dependent coupling constants have to be com- 
patible with the existence of this state. This is possi- 
ble only if ( J^j -I- -) > for any r. The conditions 
that ensure the existence of a factorized ground state are 
Jr,2 — ^"^"i; ^-iid Jr.i — lrJ\i-, and the factorizing field 
amounts to hp = (1 

J a TQ 
1 — Z^r ^r,l- 

hybrid, way to 

is to fix, 
and (Jf,2 + Jf,2) 

1=0 1^41+1/ 1^41+2/ 

ground state. This implies 

rJ^ 2 ~ i^J^ ]^ , J J. ^ — '~^rJ\ and 



^/W+jfTW+Jf). with 



The second, 
tiferromagnetism 

{Jh + Jli) > 

ing l^f) = « 
as a factorized 



introduce an- 
for example 
< 0, obtain- 

T 



j) \^4l+3 



) 1^4^+4) 



•Jr,2 



(-1) 



leading for the factorized 
with 



jx,y 



field to hp = (1 + k) VW+JFHJ^ 

In summary, we studied the zero-temperature phase 
diagram of the dimerized XY chain in a transverse field. 
We discussed the existence of a fully unentangled ground 
state, which depends on whether the parameters of the 
system satisfy given property. Furthermore, we showed 
the role of the factorizing field inside the general solu- 
tion of the model. It represents a border line between 
two separate symmetry-broken regions in the space of 
the Hamiltonian parameters. In analogy with the homo- 
geneous case, where these two regions are characterized 
by qualitatively different types of entanglement, namely 
parallel and antiparallel entanglements we expect that 
the same transition could take place also in our system. 
On the other hand, also when the Hamiltonian parame- 
ters are not compatible with ground-state factorization, 
there is a value of the field which separates the two re- 
gions. Finally, we extended the search for ground state 
factorization to more general dimerized models. 
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